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MOLECULAR CONSIDERATIONS IN THE 
TYLOPHORINAE K. SCHUM. 
(APOCYNACEAE-ASCLEPIADOIDEAE) 


S. LIEDE*, A. TAUBER* & J. ScHNEIDTT 


More than 30 collections of Tylophora from all over its geographic range were analysed 
together with representatives of all other genera in subtribe Tylophorinae except 
Rhyncharrhena. Regions of the chloroplast genome (trn T-L spacer, trnL intron and 
trnL-F spacer) as well as the nuclear genome (internal transcribed spacer including 5.8S 
of nuclear ribosomal DNA) were sequenced. Both the cpDNA dataset and the rDNA 
dataset yield a number of stable clades, but a completely resolved and stable phylogeny 
Is not obtained from either alone or by combining them. The small East African and 
Arabian genera Blyttia, Diplostigma, Goydera and Pleurostelma form a well-supported 
clade in both datasets. Pentatropis is always recognized as monophyletic. Vincetoxicum is 
not clearly separable from Tylophora, and probably covers the temperate species, while 
Tylophora covers those in the tropics and subtropics. 


Keywords. Apocynaceae, Asclepiadaceae, Asclepiadoideae, ITS, molecular systematics, 
trnT-L spacer, trnL intron, trnL-F spacer, Tylophorinae. 


INTRODUCTION 


The position of Tylophora R. Br., a genus of the Old World tropics comprising 
100—150 species, has long been misunderstood because the tribe-defining character 
of pollinia orientation (Bruyns & Forster, 1991; Liede, 1996a) is expressed somewhat 
ambiguously with usually horizontal, rarely slightly erect or hanging pollinia. 
Swarupanandan et al. (1996) showed clearly that the morphology of the stamens, 
however, conforms to that of the Asclepiadeae. The flowers of Tylophora and its 
relatives are usually rather small and simple for members of the Asclepiadoideae; the 
corona is in most cases purely or at least predominantly (e.g. T. anomala) staminal 
and attached to the back of the anthers. This simplicity has probably misled research- 
ers into constructing a link to the Marsdenieae or at least to consider Tylophora as 
a primitive genus in the Asclepiadeae. (For a thorough discussion of morphological 
characters see Omlor, 1998.) However, the position of Tylophora and Vincetoxicum 
in the matK (Civeyrel et al., 1998) and rbcL (Sennblad, 1997) trees points to a 
derived position for the Tylophorinae in the Asclepiadoideae-Asclepiadeae and thus 
supports the interpretation that floral structure is of secondary simplicity. 

The position of Tylophora as sister group to the speciose temperate genus 
Vincetoxicum Wolf has been deduced from chemical characters (Liede, 1996b) and 
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confirmed by molecular analysis of rbcL (Sennblad, 1997) and matK (Civeyrel et al., 
1998). Liede (2001) has shown that the Tylophorinae (Biondia Schltr., Blyttia Arn., 
Diplostigma K. Schum., Goydera Liede, Pentatropis R. Br., Pleurostelma Baill., 
Tylophora, Vincetoxicum and probably Rhyncharrhena F. Muell.) constitute a mono- 
phyletic, very closely interrelated subtribe, only very distantly related to the 
Astephaninae s. str., with which it was previously combined for morphological 
reasons, namely long, slender, non-verrucose hairs on the adaxial corolla surface 
and clear sparse latex (Liede, 1994, 1997). However, Forster (1992) included some 
species in Tylophora that contain milky latex, and reports on latex colour for some 
other taxa in the Tylophorinae (e.g. Vincetoxicum carnosum; Forster, 1988) are 
contradictory. 
This paper aims to answer the following questions: 


_ 


Is the present circumscription of Tylophora correct or should Tylophora or some 
of the small Tylophorinae genera (Blyttia, Diplostigma, Goydera, Pentatropis and 
Pleurostelma) be included in Vincetoxicum? 


N 


Is there a recognizable subgeneric structure in Tylophora? 


Lu 


Have the species with milky latex rightfully been included in Ty/ophora (and the 
Tylophorinae) and if so, are they related or does this trait represent parallel evol- 
utionary events? 


MATERIALS AND METHODS 
Taxa 


Thirty-two accessions in the Tylophorinae sensu Liede (2001) were analysed, among 
them 21 species of Tylophora from Africa, Asia and New Caledonia. Outgroup choice 
was based on the results of Liede (2001); a member of the Asclepiadinae 
(Gomphocarpus R. Br.) and two of the Old World Metastelminae (Cynanchum L. and 
Schizostephanus Hochst. ex Benth. & Hook.) were selected (see Table 1 for species 
names and authors). 


DNA extraction and PCR 


DNA was isolated from fresh or dried leaf tissue according to the method of Doyle 
& Doyle (1987). PCR primers and protocols for the chloroplast DNA (cpDNA) 
trnT-L and trnL-F spacers and the trnL intron correspond to Taberlet er al. (1991). 
The entire internal transcribed spacer (ITS) region including 5.8S of ribosomal 
DNA (rDNA) was amplified using the flanking primers ITS4 and ITS5 following a 
modified protocol based on Baldwin (1992) described in Meve & Liede (2001). 


Data analysis 


cpDNA sequences were pre-aligned with Perkin Elmer Sequence Navigator Version 
1.0.1; the alignment was then adjusted manually. The alignment comprised 35 taxa 
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TABLE 1. Voucher and locality information for plant material used in this study 
EMBL 
Accession No. EMBL 
trnT-L spacer Accession 
trnL intron No. ITS 
Species Origin Voucher trnL-F spacer (incl. 5.8S) 
Outgroups 
Gomphocarpus South Africa Nicholas 2829 (UDW) AJ290875 
physocarpus AJ290876 AJ320446 
E. Mey. AJ290877 
(Asclepiadeae- 
Asclepiadinae) 
Cynanchum ellipticum South Africa Liede 2933 (UBT) AJ290845 
(Harv.) R.A. Dyer AJ290846 AJ320444 
(Asclepiadeae- AJ290847 
Metastelmatinae) 
Schizostephanus alatus Kenya Noltee s.n. sub IPPS AJ410247 
Hochst. ex 8111 (UBT) AJ410248 AJ320451 
K. Schum. AJ410249 
(Asclepiadeae- 
Metastelmatinae) 
Ingroup (Asclepiadeae-Tylophorinae sensu Liede, 2001) 
* Biondia henryi (Warb. China Deng 90203 (MO) AJ410190 
ex Schltr. & Diels) AJ410191 - 
Tsiang & P.T. Li AJ410192 
Blyttia fruticulosum Kenya Liede & Newton 2946 AJ410193 
(Decne.) D.V. Field (UBT) AJ410194 AJ320443 
AJ410195 
Diplostigma canescens Kenya Liede & Newton 3214 AJ410199 AJ320445 
K. Schum. (UBT) AJ410200 
AJ410201 
Goydera somaliense Somalia Thulin & Bashir 6882 AJ410208 AJ320447 
Liede (UPS) AJ410209 
AJ410210 
Pentatropis Madagascar Liede 2749 (UBT) AJ410235 
madagascariensis AJ410236 AJ320448 
Decne. AJ410237 
Pentatropis nivalis Kenya Meve 949 (B, MSUN, AJ410238 
(J.F. Gmel.) UBT) AJ410239 AJ320449 
D.V. Field AJ410240 
& J.R.I. Wood 
Pleurostelma cernuum Tanzania Liede & Meve 3377 AJ410241 
(Decne.) Bullock (B, UBT) AJ410242 AJ320450 
AJ410243 
Tylophora anomala Cameroon Meve 916 (K, MSUN) AJ410250 
N.E. Br. AJ410251 AJ320452 
AJ410252 
Tylophora apiculata Kenya Robertson 7016 (UBT) AJ410253 
K. Schum. s.l. AJ410254 AJ320453 


AJ410255 
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TABLE 1. (Cont'd) 


EMBL 

Accession No. EMBL 
trnT-L spacer Accession 
trnL intron No. ITS 


Species Origin Voucher trnL-F spacer (incl. 5.8S) 
Tylophora biglandulosa New Veillon 7899 (P) AJ320401 AJ320454 
Endl. Caledonia AJ320402 
AJ320403 
Tylophora conspicua Tanzania Liede & Meve 3366 AJ320404 
N.E. Br. (B, UBT) AJ320405 AJ320455 
AJ320406 
Tylophora coriacea Mauritius Bernardi s.n. (BR, AJ320407 AJ320456 
(Decne.) Marais UBT) AJ320408 
AJ320409 
Tylophora flanaganii South Africa Nicholas 2839 (UDW) AJ410256 
Schltr. AJ410257 AJ320457 
AJ410258 
*Tylophora flexuosa Borneo Schneidt JS95-84 AJ320410 
R. Br. (ABD) AJ320411 E 
AJ320412 
Tylophora flexuosa Philippines — Schneidt JS96-40 AJ320413 
R. Br. (UBT) AJ320414 AJ320458 
AJ320415 
Tylophora perrottetiana Philippines — Schneidt JS96-31 AJ320416 
Decne. (UBT) AJ320417 AJ320459 
AJ320418 
Tylophora perrottetiana Philippines Liede 3252 (UBT) AJ290915 
Decne. AJ290916 AJ320460 
AJ190917 
Tylophora heterophylla Kenya Liede & Newton 3155 AJ410259 
A. Rich. (UBT) AJ410260 AJ320461 
AJ410261 
Tylophora hirsuta India Chaturvedi s.n. (UBT) AJ320419 
Wight AJ320420 AJ320462 
AJ320421 
Tylophora indica India Bruyns s.n. (MSUN) AJ410262 
(Burm.f.) Merrill AJ410263 AJ320463 
AJ410264 
Tylophora oblonga Cameroon Meve 915 (B, MSUN, AJ320422 AJ320464 
N.E. Br. UBT) AJ320423 
AJ320424 
Tylophora parviflora Philippines Liede 3290 (B, UBT) AJ320425 
(Wight) Meve, AJ320426 AJ320465 
Omlor & Liede AJ320427 
Tylophora sylvatica Africa (ex Valck s.n. (UBT) AJ410265 
Decne. hort.) AJ410266 AJ320466 


AJ410267 
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TABLE 1. (Cont'd) 


EMBL 

Accession No. EMBL 
trnT-L spacer Accession 
trnL intron No. ITS 


Species Origin Voucher trnL-F spacer (incl. 5.8S) 
Tylophora Kenya Liede & Newton 3200 AJ320428 
tenuipedunculata (MSUN) AJ320429 AJ320467 
K. Schum. AJ320430 
Tylophora tenuis Borneo Schneidt JS95-110 AJ320431 
Blume (ABD) AJ320432 AJ320468 
AJ320433 
Tylophora villosa Borneo Schneidt JS95-39 AJ320434 
Blume (ABD) AJ320435 AJ320469 
AJ320436 
Tylophora villosa Philippines — Schneidt JS96-109 AJ320437 
Blume (ABD) AJ320438 AJ320470 
AJ320439 
Tylophora yunnanensis China Schneidt JS96-134 AJ320440 
Schltr. (ABD) AJ320441 AJ320471 
AJ320442 
Vincetoxicum atratum China Schneidt JS96-137 AJ410268 
Morr. & Decne. (ABD) AJ410269 AJ320472 
AJ410270 
Vincetoxicum carnosum Borneo Schneidt JS95-97 AJ410271 
Benth. (ABD, L) AJ410272 AJ320473 
AJ410273 
Vincetoxicum Germany Meve 970 (UBT) AJ410274 
hirundinaria Medic. AJ410275 AJ320474 
AJ410276 
Vincetoxicum stocksii Pakistan Ali & Khatoon s.n. AJ410277 
S.I. Ali & (GA) AJ410278 AJ320475 
S. Khatoon AJ410279 


*ITS sequence not available. 


and 1805 characters; 80 data cells were unknown, affecting the beginning of the 
trnT-L spacer in Biondia henryi and T. apiculata, the end of the trnT-L spacer in T. 
yunnanensis, the beginning of the trnL intron in B. henryi, the end of the trnL intron 
in P. nivalis, and the end of the trnL-F spacer in T. coriacea and T. hirsuta. 

Indels were coded as 'missing characters; parsimony-informative indels were 
coded separately as present or absent following the ‘simple gap coding’ method of 
Simmons & Ochoterena (2000). Indels resulting from a different number of the same 
base in a chain of four or more repeats of the same base were not coded separately, 
as the length of these chains has been found to vary even within the same species 
(Liede, unpubl. data). This affected bp 565—570 and 590—601 in the trn T-L spacer, bp 
295—320 in the trnL intron, as well as two of the three potentially informative indels 
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in the trnL-F spacer. The 11 separately coded indels were added to the matrix in a 
second analysis. The alignment is available from the senior author and can be viewed 
in TreeBase (SN940; Sanderson et al., 1994). 

Not all rDNA sequences could be aligned unambiguously (see Appendix; TreeBase 
SN941). Manual alignment of sequences produced a lower number of parsimony- 
informative characters than computer programs tested for different algorithms and 
settings (Jotun-Hein algorithm in ‘Megalign’ of DNAStar (Lasergene) and Clustal 
Algorithm in Perkin Elmer Sequence Navigator Version 1.0.1) and was thus con- 
sidered the most conservative approach. Even so, both ITS! and ITS2 contain large 
areas of more or less ambiguous alignment. Excluding those areas introduces a 
subjective view as to which areas are sufficiently well aligned to include and which 
ones need to be excluded. Nevertheless, in a second analysis, the 135 most ambigu- 
ously aligned characters were excluded, 68 in ITS1 and 67 in ITS2 (see Appendix). 
This selective approach was preferred over a general exclusion of missing or ambigu- 
ous characters, which would have reduced the number of parsimony-informative 
characters to 73. Separate coding of indels was not carried out because of the hetero- 
geneity of the dataset. 

As a partition homogeneity test did not show significant discordance between the 
cpDNA and the rDNA dataset, a combined analysis was performed. Again, all 
sequence characters were analysed, and, finally, the 11 separately coded indels for 
the cpDNA dataset were added. The two taxa for which no rDNA data were avail- 
able (Biondia henryi and Tylophora flexuosa JS95-84) were coded unknown for all 
rDNA data and included in the combined dataset. 

Sequence analysis, phylogenetic analysis and tests for clade support were per- 
formed using PAUP version 4.0d65a (PPC) (Swofford, 1998) on a Macintosh G3 
Powerbook. 

Sequence divergence among taxa was calculated using the 'show pairwise distance 
option for the trnT-L spacer, the trnL intron and the trnL-F spacer, as well as for 5.8S, 
ITSI and ITS2 (excluding the end of 18S and the beginning of 26S). The beginning 
and end of ITS! and ITS2 were determined in comparison with Asclepiadoideae 
sequences submitted to EMBL (AJ402152-AJ402162, Meve & Liede, 2001). 

For parsimony analysis, heuristic searches were conducted on all datasets (addition 
sequence random, 100 replicates, TBR branch swapping, ‘MulTrees’ on, ‘steepest 
descent’ off ). Bootstrap search (1000 replicates) was conducted under the ‘full heuris- 
tic’ search option, swapping algorithm was set to ‘TBR’, ‘MulTrees’ on, ‘steepest 
descent’ off, and ‘maxtrees’ set to 1000. Jackknife resampling (1000 replicates) was 
set to 50% deletion, and ‘Jac’ resampling; the other settings were identical to the 
bootstrap settings. 

Finally, a maximum likelihood analysis was performed on both single and com- 
bined datasets under the HK Y85 model (Hasegawa et al., 1985), using the standard 
settings as implemented in PAUP 4.0b5 (PPC) (Swofford, 1998). As computational 
limits were a serious concern, all taxa representing duplicate accessions of the same 
taxon and closely related species were excluded, resulting in 31 taxa for the cpDNA 
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dataset and 30 for the rDNA and the combined dataset. Ten replicate searches were 
conducted for each dataset with the heuristic search algorithm as implemented in 
PAUP 4.0b5 with random addition sequence and TBR options invoked. 


RESULTS 


Sequence characteristics for both datasets are summarized in Table 2. 

Sequence divergence is higher in ITSI, with 1.5—-15.7% between ingroup taxa and 
6.2-17.3 between ingroup and outgroup taxa, than in ITS2, with 0.8-13.0% and 
3.0—15.196, respectively. For the cpDNA, the ¢rnL intron has the lowest sequence 
divergence, with 0.0-1.3% between ingroup taxa and 0.4-1.7% between outgroup 
and ingroup taxa, while the trnL-F spacer has the highest sequence divergence, with 
0.0-3.0% and 0.5—3.396, respectively. Thus, sequence divergences in rDNA are 
between five and ten times higher than in cpDNA. 

As an analysis of the sequence characters of the cpDNA dataset (independent of 
separately coded indels) resulted in more than 60,000 most parsimonious trees and 
thus exceeded computer memory capacity, the topology of the strict consensus trees 
was obtained excluding uninformative characters. Both the strict consensus tree 
(Fig. 1A) and the ‘semistrict’ consensus tree show a number of well-supported clades, 
but no resolution between them. The addition of the indels hardly changes these 
clades, only adding 7: parviflora to the Pentatropis clade, but without any support 
(Fig. 1A). 

The maximum likelihood analysis under the HK Y85 model (Hasegawa et al., 
1985) yielded one most likely tree (-Ln likelihood =3687.90, Fig. 1B). The 
Tylophorinae are clearly monophyletic, but comprise several unresolved groups: the 
small African genera, Pentatropis—T. parviflora, T. sylvatica-T. apiculata-T. tenuipe- 
dunculata, and the remainder of Tylophora—Vincetoxicum with subgroups corre- 
sponding largely to the clades of the parsimony analysis. 


TABLE 2. Sequence characteristics 


trnT-L trnL trnL-F 
spacer intron spacer ITSI ITS2 5.8S 
Aligned total length (bp) 877 533 384 322 291 158 
Length range (bp) 790-832 504—533 354-374 194-242 217-247 158 
Length mean (bp) 800.80 517.40 369.14 225.27 234.36 158 
Number of parsimony- 7 3 1 _ _ _ 
informative indels* 
Sequence divergence 0.0-1.9 0.0-1.3 0.0-3.0 (0)t 1.5- (0) 0.8- 0.0-3.8 
(ingroup) (%) 15:7 13.0 
Sequence divergence 1.0-2.9 0.4-1.7 0.5-3.3 6.2-17.3 3.0-15.1 0.0-2.5 
(ingroup/outgroup) 
(%) 


*Not considering variable lengths of chains >4 bp. 
TBetween T. perrottetiana SL3252 and JS96-40 (see text). 
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Fic. 1. cpDNA dataset. A, Strict consensus of more than 60,000 most parsimonious trees in 
one island (17 trees without uninformative characters) resulting from parsimony analysis of 
all sequence characters and strict consensus of more than 60,000 most parsimonious trees in 
one island (79 without uninformative characters) resulting from parsimony analysis of all 
sequence characters and indels. Numbers above branches indicate bootstrap/jackknife percent- 
ages for the analysis without indels, those below the branches for the analysis with indels. 
Dotted lines indicate clades not retrieved without indels. -/- denotes values below 50%. B, 
Maximum likelihood tree of cpDNA sequence characters (Ln = 3687.9). 


For the rDNA dataset, difficulties in alignment are reflected in low indices for 
measuring homoplasy (Table 3). However, the clades retrieved are identical to those 
retrieved from the cpDNA dataset, only the most basal taxa of single clades (T. 
yunnanensis, T. heterophylla) are unresolved for the rDNA (Fig. 2A). Not even the 
Cynanchum-Schizostephanus clade is present in all most parsimonious trees, and 
therefore it does not occur in the strict consensus, despite bootstrap and jackknife 
values of more than 80% (Fig. 2A). 

Excluding the ambiguously aligned characters of the rDNA dataset retrieves all 
clades found with all rDNA characters, plus the Cynanchum-Schizostephanus clade, 
the T. conspicua-T. oblonga clade and the V. hirundinaria—V. stocksii clade (Fig. 2B). 
The most obvious difference is that the Cynanchum-Schizostephanus clade is now a 
member of the ingroup, with T. tenuis in a basal position; this unexpected position 
for T. tenuis is without bootstrap or jackknife support. If the monophyly of the 
Tylophorinae is enforced by constraint, the most parsimonious trees are three steps 
longer (/ — 411) than without constraint, a result considered not significant in a 
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= — — — — — — Biyttia friticui 19176 Blyttia fruliculosum 
28 See Diplostigma canescens 
98/96 E Diplostigma canescens 95/92 
Goydera somaliense Goydera somaliense 
Tylophora biglandulosa Tylaphora biglandufosa 
TT unu Tylophora anomale 
sso [— — — — — —— lylophora anomala 86/76 : 
97195 LL Tylophora coriacea 97/94 Tylophora coriacea i 
Tylophora flanaganii Tylophora flanaganii 
"S [OO Tylophora Mais mm Fs Tylophora ax 
071 nss e Tyiophora indica — — ] aia Tylophora indica, 
— 93691 Tylophora villosa JS95 39 95/88 Tylophora villosa J595 39 
A Tylophora villosa JS96 109 B Tylophora villosa JS96 109 


Fic. 2. rDNA dataset. A, Strict consensus of the 12 most parsimonious trees in three islands 
resulting from parsimony analysis of all sequence data. The dotted line indicates a clade not 
present in the strict consensus, but with a high bootstrap/jackknife support. B, Strict consensus 
of the 110 most parsimonious trees in one island resulting from parsimony analysis of sequence 
data without 135 ambiguously aligned characters. Numbers indicate bootstrap/jackknife 
percentages. 


PTP (Permutation Test for the Existence of Phylogenetic Structure), but significant 
in a T-PTP (Topology Dependent PTP) as implemented in PAUP. If ITSI is analysed 
alone, with ambiguous characters removed (228 characters, 48 parsimony- 
informative ones), all ingroup clades and Cynanchum-Schizostephanus are unre- 
solved; if ITS2 is analysed alone with the ambiguous characters removed (250 charac- 
ters, 54 parsimony-informative ones), all ingroup clades and  Cynanchum- 
Schizostephanus remain unresolved, but T. tenuis retains a basal position. 

The maximum likelihood trees of the rDNA data, with and without ambiguous 
characters (— Ln likelihood — 3292.67 and 3289.86, respectively), show no additional 
information; in the tree without ambiguous characters, 7. tenuis again takes a basal 
position. 

The combined dataset was run with all rDNA characters and with the ambiguous 
rDNA characters excluded, both with and without addition of the separately coded 
cpDNA indels (Table 3). In the analysis with all rDNA characters, but without 
indels, Pentatropis is basal in the Tylophorinae, while all other Tylophorinae are 
unresolved, forming largely the same subclades as in the previous analyses, though 
with higher bootstrap/jackknife support (Fig.3A). Adding the indels makes 
Pentatropis and T. parviflora sister to T. sylvatica-T. apiculata-T. tenuipedunculata, 
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Gomphocarpus physocarpus Gomphocarpus physocapus 
Cynanchum ellipticum Cynanchum ellipticum 
90/90 L Schizostephanus alatus 90/87 Schizostephanus alatus 
Pentatropis madagascariensis [-————~ Pleurostelma cernuum 
100 d Wale Wës 94/9 Biytia iruticulosum 
——  -NTN——---.T Biontie henryi 79179 r— Diplostigma canescens 
Tylophora heterophylla 100/99. Goydera somaliense 
Tylophora parviflora + Tylophora parviflora 
Tylophora tenuis + Pentalropis madagascariensis 
ee atratum z mE Pentatropis nivalis 
‘incefoxicum camosum KÉ Tylophora sylvatica 
86/88 [——— Tylophora conspicua 7878 ;—— Tylophora apiculata 
82/60 — Tylophora oblonga 85/87 Tylophora tenuipedunculata 
Tylophora sylvatica E Biondia henryi 
7507 r— Tylophora apiculata 8292] ————F pÑP— ——— "ylophora helerophylla 
82/83 —— Tylophora tenuipedunculata [ ka,  —- Ywiceloxicum atratum 
Pleurostetma cemuum F — , a ə yp =Á “ — —— Vincetoxicum camosum 
= $382 Biyttia fruliculosum Tylophora conspicua. 
TIG [—— Uiplostigma canescens 60/58 es Tylophora ablonga 
100/100——— Goydera somaliense Vincetoxicum hirundinana 
Tylophora biglandulosa 58/54 —— Vincetoxicum stocksii 
59/61 SE anomala -1* Tylophora yunnanensis 
70174 lophora coriacea E Tylophora hirsuta 
T5IT5 Tylophora flanaganii Tylophora indica 
Tylophora flexuosa JS95 84 Mn Tylophora vilose J885 39 
Tylophora flexuosa JS96 40 SR 88/85 Tylophora villosa JS96 108 
64/64. Er Tylophora perroftetiana JS96 31 ——TWTYÑTylophora biglandulosa 
Tyfophora perrottetiana SL3252 kaal Tylophora anomata 
rn Vincetoxicum hirundinaria TES r Tylophora coriacea 
57/52 Vincetoxicum stocksii i = T6i72 — Tylophora flanaganii 
F3 r~ — Tylophora yunnanansis r Tylophora tenuis 
E Tylophora hirsuta T Tylophora flexuosa JS95 84 
98551 Tytophora indica Tylophora flexuosa JS96 40 
A 91/82 Tylophora villosa JS95 39 B 64164 r—— lWophora perrottetiana JS96 31 
90/88 Tylophora villosa JS96 109 89/88 —— Tylophora perroitetiana SL3252 





Fic. 3. Combined cpDNA and rDNA dataset, all rDNA data included. A, Strict consensus of 
the 396 most parsimonious trees in one island resulting from parsimony analysis of all sequence 
data. -/- denotes values below 50%. B, Strict consensus of the 33 most parsimonious trees in 
one island resulting from parsimony analysis of all sequence data and the indels. -/- denotes 
values below 50%. 


and the two together to the small African genera (Fig. 3B). This whole group is 
sister to the remaining Tylophorinae, in which the T. flexuosa clade is sister to the 
T. biglandulosa-T. anomala-T. coriacea-T. flanaganii clade. All this additional struc- 
ture, however, is without any support. 

Excluding the ambiguously aligned rDNA characters, Pentatropis and the T. 
sylvatica-T. apiculata-T. tenuipedunculata group are basal to the remaining 
Tylophorinae if the indels are not added, otherwise the same clades as in the single 
analyses are retrieved (Fig. 4A). 

Maximum likelihood analysis of the combined dataset without ambiguously 
aligned rDNA characters places the T. sylvatica-T. apiculata-T. tenuipedunculata 
group basal to all Tylophorinae, and leaves a Pentatropis-T. conspicua-T. oblonga 
group unresolved with a 7. heterophylla-T. parviflora-small African genera group 
and a group comprising the remainder of Tylophorinae in three subgroups (Fig. 4B). 


DISCUSSION 
Clade structures 


In the small genera, the outcome of Pentatropis as a distinct clade is unambiguous, 
and it occupies a rather isolated position in the ingroup. A basal position for 
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Gomphocarpus physocarpus r——— PPO — — Gomphocatpus physocarpus. 
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Tylophora yunnanensis Ty 
'yfophora flexuosa JS96 40 
62/57 Vincetoxicum hirundinaria H i 
+ ^ Tylophora perrottetiana SL3252 
ators 61/57 Vincetoxicum stocksit " 
Tylophora hirsuta Tylophore biglandulosa 
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Fic. 4. Combined cpDNA and rDNA dataset, ambiguously aligned rDNA characters excluded. 
A, Strict consensus of the 5232 most parsimonious trees in two islands resulting from parsi- 
mony analysis of all sequence data and of 444 trees in two islands resulting from parsimony 
analysis of all sequence data and the indels. Numbers above branches indicate bootstrap/ 
jackknife percentages for the analysis without indels, those below the branches for the analysis 
with indels. Dotted lines indicate clades not found with indels. -/- denotes values below 50%. 
B, Maximum likelihood tree resulting from analysis of cpDNA and rDNA, ambiguously 
aligned rDNA characters excluded (Ln — 7465.55). 


Pentatropis in the Tylophorinae, however, is suggested only by all rDNA data 
(Fig. 3A). Vincetoxicum carnosum, despite considerable morphological similarity 
with Pentatropis (Liede, 1994), is always very isolated in the Tylophorinae and does 
not join the Pentatropis clade. The two species of Pentatropis are very closely related 
and show only three differences in the cpDNA dataset and six in the rDNA dataset. 
Floral structure 1s also very similar, so that possibly the much more tender, small- 
leaved Madagascan P. madagascariensis might be considered a subspecies under P. 
nivalis (see Liede & Meve, 2001). 

Very constant throughout is the Diplostigma- Goydera clade, indicating a very close 
relationship between these two monotypic East African genera. The morphological 
analyses by Liede (1994) and Bruyns (1999) did not suggest this affinity for the 
morphologically very unusual Goydera (Liede, 1993). In all analyses, the four small 
African genera form a well-supported clade. While the cpDNA data suggest a 
Blyttia-Pleurostelma clade as sister to Diplostigma—Goydera (Fig. 1), the rDNA data 
and the combined datasets place Blyttia as sister group to the Diplostigma—Goydera 
clade and Pleurostelma basal to the whole clade. In both morphological analyses 
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(Liede, 1994; Bruyns, 1999), Pleurostelma, mostly because of its long stylar head, 
was placed basal to the (Oncinema)—Astephanus—Microloma clade that has been 
shown to be of different affinity (Liede, 2001). 

Of the Tylophora species analysed, three African clades can be discerned: T. sylvat- 
ica always comes out with T. tenuipedunculata and T. apiculata, though the basal 
position in the clade is sometimes occupied by 7. apiculata, and sometimes by T. 
sylvatica. Tylophora sylvatica and T. tenuipedunculata share small flower size and 
show similar corona structure, but T. tenuipedunculata has rather reduced inflor- 
escences. Tylophora apiculata shares with T. sylvatica the more floriferous inflor- 
escences, but in addition has larger flowers. The second clade of African Tylophora 
species comprises 7. anomala, T. coriacea and T. flanaganii. Morphologically, the 
well-developed, erect, firm corona lobes and the large pollinaria with hanging pollinia 
characterize this group of species. In the rDNA and the combined datasets, this 
clade is associated with the New Caledonian T. biglandulosa, a species with white 
milky latex (Forster, 1991). Forster (1991) is thus supported in his view that there 
is no reason to separate T. biglandulosa from Tylophora as Hybanthera Endl. The 
third clade of African Tylophora species, T. conspicua and T. oblonga, is retrieved 
in all analyses except for that which includes all rDNA data, but the clade is more 
strongly supported by cpDNA data than by rDNA data. Morphologically, corona 
structure is very different in the two species, but both possess unusually large flowers 
and broad, apically twisted corolla lobes. T'ylophora oblonga, originally in the mono- 
typic genus Oncostemma K. Schum. because of its unique corona structure (Meve 
& Omlor, 2000), has thus correctly been included in Tylophora. Tylophora conspicua 
possesses whitish, milky latex. 

The Asian Tylophora species fall largely into two clades. The first consists of 7: 
indica and the two specimens of T. villosa, indicating a relationship supported by 
various morphological traits (Schneidt, 1999). This clade is always basally joined by 
T. hirsuta, which has also been considered a member of this complex for morphologi- 
cal reasons. For the cpDNA dataset and both combined datasets, this clade is basally 
joined by the Chinese 7: yunnanensis and also by Vincetoxicum hirundinaria and V. 
stocksii. Only in the rDNA dataset analyses (Fig. 2) do T. yunnanensis and the V. 
hirundinaria-V. stocksii clade remain unresolved; in the analysis of all rDNA data, 
the latter two species also remain unresolved. Flower structure also supports an 
association of these species with the T. indica group. 

The second clade comprises most members of the T. flexuosa group. Schneidt 
(1999) proposes four varieties in this complex, of which three were included in the 
present analysis. The two Philippine accessions of T. perrottetiana (SL3252 and 
JS96-31) are almost identical in the regions studied, with no difference in their rDNA 
sequences and only three in the trnT-L spacer. Characteristic of this species are two 
deletions in the trnL intron, one being unique in the dataset and the other unique 
in the ingroup. The other specimen from the Philippines (JS96-40) is clearly more 
distinct, with 20 differences from JS96-31 in the rDNA dataset. Owing to insufficient 
material no rDNA sequence could be obtained from the Bornean specimen 
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(JS95-84). This differs from JS96-40 in nine sites in the cpDNA dataset. This number 
of differences is the same as between the two specimens of T. villosa, one from 
Borneo (JS95-39) and one from the Philippines (JS96-109), with five sites differing 
in both the cpDNA dataset and the rDNA dataset. The last specimen of the T. 
flexuosa complex, JS95-110, represents T. tenuis, also proposed as a variety of T. 
flexuosa by Schneidt (1999). However, this specimen joins the T. flexuosa clade only 
in the analysis combining all data (Fig. 3B) and the maximum likelihood analysis 
of all sequence data without the ambiguously aligned characters (Fig. 4B); it is 
usually unresolved or weakly attached to other species. In the analysis of rDNA 
data without ambiguously aligned characters, T. tenuis is actually basal to all 
Tylophorinae and Cynanchum—Schizostephanus. Schneidt (1999) distinguished T. per- 
rottetiana from T. flexuosa by corona differences, and from T. tenuis by differences 
in inflorescence structure. These differences in corona structure may be best recog- 
nized at varietal rank. However, differences in inflorescence structure indicate differ- 
ences in species level, as none of the molecular datasets gives any indication that T. 
tenuis is not an independent species, despite its floral similarities to T. flexuosa. 

The remaining species of Tylophora (the aberrant Philippine T. parviflora (Meve 
et al., 2002) and the African T. heterophylla) and Vincetoxicum (V. atratum) are 
unambiguous members of the Tylophorinae, but show no constant affinities to any 
clade. No rDNA sequence was available for Biondia, and the cpDNA sequences 
confirm its membership of the Tylophorinae, but give no indication as to its position 
within the subtribe. 


Evolutionary implications 


Both molecular regions recognize the Tylophorinae sensu Liede (2001) as a well- 
defined monophyletic group with high bootstrap and Jackknife support. While the 
poor resolution of the cpDNA results is easily understood because of lack of variation 
between the sequences, the problems with the rDNA dataset are unexpected. 
Sequence divergence is much lower in rDNA datasets from other Asclepiadoideae 
genera and genera complexes (e.g. < 10% in the succulent Stapelieae, Meve & Liede, 
unpubl. data; «1596 in Cynanchum s. str., Liede, unpubl. data), but not unusual for 
this region in other angiosperms (e.g. Streptocarpus, Móller & Cronk, 1997; 
Heliantheae, Urbatsch et al., 2000). However, the problem of this dataset lies in the 
high number of very variable indels interspersed with short, hardly variable regions. 
This probably reflects conservation of secondary structure, but not of primary 
sequence (Coleman et al., 1998). 

Therefore, both datasets result in groups of species clustering together, though 
basal resolution is undetectable. Under almost all settings, the 'semistrict' consensus 
does not show any additional structure compared with the strict consensus. The 
clades regularly retrieved largely coincide in both rDNA and cpDNA analyses, and 
most of them are supported in one of the single datasets or one of the combined 
ones by satisfactorily high bootstrap/jackknife values (> 80%), so that they can be 
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regarded with some confidence as monophyletic units. These monophyletic units are 
composed largely of species of close geographic distribution, though the New 
Caledonian T. biglandulosa shows more affinities to the African T. anomala-T. 
coriacea-T. flanaganii clade than to any of the Asian clades. 

In terms of an evolutionary scenario, poor differentiation of the chloroplast spacer 
regions and the trnL intron and the patchwork indels of ITS1 and ITS2 point to 
very rapid evolution of geographical groups from a common, possibly widely distrib- 
uted, ancestor. As previous studies of the Apocynaceae as a whole have shown 
(Sennblad, 1997; Civeyrel et al., 1998), Tylophora is one of the most derived genera 
of the Asclepiadoideae, so that longevity of the group as reason for high ITS variation 
is unlikely. However, Tylophora and Vincetoxicum are the only genera in 
Asclepiadoideae for which autogamy and apomixis have been demonstrated 
(Chaturvedi, 1988; Naumova, 1992; Lumer & Yost, 1995). Furthermore, their 
smooth seeds are very effectively wind dispersed, and the normally small, rather 
simple, but often very nectariferous flowers are commonly visited by small insects 
not covering long distances. Thus, effective isolation mechanisms together with good 
dispersal capacity might have led to the pattern observed. On Krakatau, three out 
of 15 Asclepiadoideae reported 111 years after the catastrophic eruption are members 
of Tylophora (Bush et al., 1995). The rapid spread of Vincetoxicum in the Americas 
(e.g. Pringle, 1973; Sheeley & Raynal, 1996) and of Tylophora indica in the West 
Indies (Schlechter, 1899) further demonstrates the vigour of this group. Possibly 
several groups evolved at different times in some regions, which would explain the 
African clades having similar distribution but not close relationships, indicating 
independent colonization events in Africa. The Asian clades also appear distantly 
interrelated. The association of the type species of Vincetoxicum (V. hirundinaria) 
with the T. indica clade suggests that Tylophora comprises the tropical and 
Vincetoxicum the temperate representatives of the same complex. This view is sup- 
ported by numerous intermediate species in the Far East, which have been shuffled 
back and forth between Vincetoxicum and Tylophora. 

The two possible approaches to translate this phylogeny into a classification both 
create more problems than they solve. First, all the taxa studied here could be 
integrated into a large genus Vincetoxicum. To lump large and complex genera such 
as Vincetoxicum, Tylophora and Biondia would create a genus relatively easy to 
delimit, but very difficult to differentiate at the species level. Small, but easily reco- 
gnizable, entities such as Blyttia, Diplostigma, Goydera and Pleurostelma would have 
to be lumped into this conglomerate as well, marking the extremes of morphological 
variation, and thus blurring the genus concept. 

The second possibility, giving each clade generic recognition, would result in 
numerous small genera of uncertain delimitation that cannot be identified morpho- 
logically. However, the four morphologically very distinct East African genera could 
then remain distinct. 

The comparatively low morphological differentiation in the whole complex and 
the slightly smaller number of name changes involved support the first possibility. 
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However, before such far-reaching changes are made, a thorough study at species 
level needs to be carried out for the whole complex. Especially in mainland Asia 
and Japan, the relationships of Biondia, Tylophora and Vincetoxicum need to be 
re-examined morphologically, karyologically and with a different set of molecular 
regions, based on the results presented here. The insufficiently known monotypic 
genera Merrillanthus Chun & Tsiang and Pentastelma Tsiang & P.T. Li possibly 
represent the extremes of morphological variation in the complex in China and 
should be included in such a study. 

Finally, the present study answers the question of whether latex colour is a reliable 
synapomorphy for the Tylophorinae. While it is certainly indicative, it 1s neither 
unique (also being found in the Ceropegieae and the Astephaninae s. str.), nor present 
throughout. Tylophora parviflora with milky, yellowish latex is always placed in the 
main clade, but shows no clear affinities to any of the constant clades. The other 
two species reported to have white milky latex (T. conspicua and T. biglandulosa) 
show no close affinity to each other nor to T. parviflora. For V. carnosum, contradic- 
tory statements have been found on herbarium labels, so that seasonal changes or 
changes between populations cannot be excluded. From these results it is likely that 
a reversal to milky latex has occurred more than once in the Tylophorinae. 
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